function qs = Elnqmu_losses(mix)

% E[ln q(mu)]

T = mix.varposterior.T;
[K D] = size(mix.centres);

Dlog2pi = D*log(2*pi);

qs = zeros(1, K);
for k = 1:K
    qs(k) = - 0.5*D - 0.5*Dlog2pi - 0.5*logdet(T(:, :, k));
end
qs = sum(qs) - qs;